Quantum measurement of hyperfine interaction in nitrogen- vacancy center 
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We propose an efficient quantum measurement protocol for the hyperfine interaction between the electron 
spin and the 15 N nuclear spin of a diamond nitrogen-vacancy center. In this protocol, a sequence of quantum 
operations of successively increasing duration is utilized to estimate the hyperfine interaction with successively 
higher precision approaching the quantum metrology limit. This protocol does not need the preparation of 
the nuclear spin state. In the presence of realistic operation errors and electron spin decoherence, the overall 
precision of our protocol still surpasses the standard quantum limit. 



I. INTRODUCTION 

The negatively charged nitrogen-vacancy (NV) center in di- 
amond is a promising solid state system for quantum compu- 
tation. The electron spin in the optical ground state of the NV 
center exhibits exceptionally long coherence time (> 350 jus) 
at room temperature. 1 This feature allows coherent manipula- 
tion and reliable readout of the state of the electron spin and 
the neighboring nuclear spins 2 - 3 in the NV center, a key tech- 
nique of diamond-based quantum computation^^ 5 - In these 
operations, the hyperfine interaction between the electron spin 
and the neighboring nitrogen nuclear spin plays an important 
role. To minimize the operation errors, an accurate estimate 
of the hyperfine interaction is desirable. 

In addition to quantum computation, the NV center is also 
a candidate for the application of quantum parameter estima- 
tion (also known as quantum metrology). Quantum metrol- 
ogy seeks quantum measurement protocols to estimate phys- 
ical parameters up to a given precision defined as 1/A 2 (with 
A being the standard deviation) using the least amount R of 
resources, which include the number of measurements, the 
total duration of the measurements, and the number of par- 
ticles involved in the measurements. The classical protocol 
utilizes the number R of repeated measurements as a resource 
and, according to the central limit theorem, gives the classical 
limit (also known as standard quantum limit or SQL) Asql = 
0(1 1 V/?). Quantum metrology aims to surpass the SQL and, 
more ambitiously, reach the quantum metrology limit (QML) 
Aqml = 0(1 1 R), the upper precision bound 1/Aq ML = 0(R 2 ) 
set by quantum mechanics. The most popular quantum mea- 
surement technique is interferometry, in which the parameter 
to be measured is recorded as a phase in the coherence of the 
system. 6-9 The exceptionally long coherence time of the NV 
center electron spin diminishes the detrimental effect of deco- 
herence on such measurements and makes the NV center an 
ideal system for quantum metrology. 10 Up to date, most of the 
measurement protocols utilize pure quantum states and sur- 
pass the SQL by creating quantum entanglement in the sys- 
tem. However, the thermal equilibrium state of the nuclear 
spins is highly mixed at room temperature. To estimate reli- 
ably the hyperfine interaction in the NV center by a pure-state 
protocol, the nuclear spins must be prepared repeatedly into a 
given pure state. Further, the number of spins as the resources 
of entanglement in a single NV center is finite, 5 so the advan- 
tage of quantum entanglement to parameter estimation is also 



limited. 

Recently, Boixo and Somma 1 1 proposed a model of mixed- 
state quantum metrology by combining the mixed-state quan- 
tum computation (also known as deterministic quantum com- 
putation with one quantum bit 12 or DQC1) with the adaptive 
Bayesian inference. This DQC1 model utilizes the total du- 
ration T (instead of large-scale entanglement 10 ) of the esti- 
mation process as a resource to approach the QML Aqml = 
0(1 /T) without creating any entanglement.^!! However, its 
application to estimate the hyperfine interaction in the NV 
center requires including the effects of noise and unintended 
dynamics. 

In this paper, we construct an efficient quantum measure- 
ment protocol to estimate the hyperfine interaction between 
the electron spin and the 15 N nuclear spin in the NV cen- 
ter. This protocol is essentially a combination of the DQC1 
model 11 and the spin-echo technique, 15 which decouples the 
dynamics driven by the hyperfine interaction from the noise 
and unintended dynamics. It does not need the preparation 
of the nuclear spin state and approaches the QML Aqml = 
0(1 /T) in the ideal case. By including realistic errors (such as 
the nuclear spin rotation error and the electron spin decoher- 
ence) in our analysis, we show that our protocol still surpasses 
the SQL under typical experimental conditions. 

The rest of this paper is organized as follows. In Sec. HI] we 
review the DQC1 model for parameter estimation and identify 
the problems in applying this model to estimate the hyperfine 
interaction in the NV center. In Sec. [Ill] we give a solution 
to these problems by combining the DQC1 model with the 
spin-echo technique. In Sec. IIVI we introduce our quantum 
measurement protocol. Sec. [V]gives the conclusion. 



II. DQC1 PARAMETER ESTIMATION IN NV CENTER 



We first review the two-qubit version of the DQC1 pa- 
rameter estimation model proposed by Boixo and Somma 11 
(Sec. Ill Al l and then identify the robustness problems arising 
from applying this model to estimate the hyperfine interaction 
in the NV center (Sec. lHBI i. 
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FIG. 1. DQC1 parameter estimation with one control qubit in the 
pure state |+> = (|0> + |1»/ V2 and one target qubit in the state p tar . 



A. Two-qubit DQC1 parameter estimation 

The two-qubit DQC1 model consists of a control qubit 
(with states {|0) , [1)}) and a target qubit (with states {|T> , ID))- 
The initial state pdqci = 1+) (+1 ® Ptar is the direct product of 
the pure state |+) = (|0) + 1 1 ))/ a/2 of the control qubit and the 
unpolarized state p tar = (|t> <TI + II) UD/2 of the target qubit, 
as shown in Fig. Q] The three Pauli operators of the control 
qubit and of the target qubit are denoted by 

*=|1><0| + |0><1|, 

?si(|i><oi-|o><iD, 

Z = |0><0|-|1><1|, 

and [fr x , fry, fr z }, respectively. The two qubits are coupled by 
the interaction 



#DQC1 = |l)<l|®ftx z . 



(1) 



This interaction makes the splitting energy u> c of the control 
qubit dependent on the state of the target qubit: <j Cj f = 9 for 
the target qubit in the spin-up state |f) and w c j = —9 for the 
target qubit in the spin-down state |X). The DQC1 parameter 
estimation 11 aims to estimate the interaction strength 9 with 
the standard deviation Ag = 0(1 /T) approaching the QML, 
where T is the total duration of the estimation process. The 
procedures are simple: the application of the two-qubit inter- 
action Hoqci for a duration r, followed by a measurement of 
X: 

• If the target qubit is in the spin-up state |f), then Hdqci 
drives the precession of the control qubit with angular 
frequency <o C) f, 



10) + ID 
V2 



'IT) 



10) + « 



ID 



V2 



®IT). 



Before the measurement, the interaction strength 9 is 
encoded as a phase e~ la> ^ T of the control qubit. The 
repeated measurements of X estimate the average value 
(X)i = cos(<y Cj |T) = cos(#r), which yields the phase. 

If the target qubit is in the spin-down state then 
^dqci drives the precession of the control qubit with 
angular frequency cj c ^, 



10) + U) 
V2 



>I4) 



|0) + e -^ T |l) 



•II). 



Before the measurement, the interaction strength 9 is 
encoded as a phase e ~ m ' iT of the control qubit. The 
repeated measurements of X estimate the average value 
(X)i = cos(w Ci |t) = cos(#r), which extracts the phase. 

Now the target qubit is in the unpolarized state, i.e., 
an equal, incoherent mixture of |f) and ID- Then 
the repeated measurements of X estimates the equally 
weighted average of and 



<*> = 5«*>r + <*h) 



cos(#r). 



A distinctive feature of the above parameter estimation 
process is the absence of any two-qubit entanglement. 14 

For a given standard deviation Ax (<k 1 under typical situ- 
ations) in estimating (X), the DQC1 model gives an estimate 
to the interaction strength 9 with a standard deviation 



(2) 



\d(X)/dd\ t\ sin(0r)| 



By regarding the duration r of the estimation as a resource, the 
QML scaling Ag = 0(1 /r) is achieved if r could be chosen 
such that | sin(#r)| w 1. However, due to the limited prior 
knowledge about 9 (the parameter to be estimated), we cannot 
always ensure | sin(0r)| m 1, especially when a small standard 
deviation Ag — > (corresponding to large t — > oo) is required. 

To address this issue, Boixo and Somma 11 quantified the 
prior knowledge about 9 by a standard deviation Aq and uti- 
lized the adaptive Bayesian inference to reduce the standard 
deviation successively. The essential idea of this approach 
can be understood qualitatively as follows. In order to ensure 
| sin(6 l r)| ss 1 and hence the QML, the largest r is roughly 
1/A . Under this restriction, the minimal standard deviation 
for the estimation of 9 is given by Eq. (f2j) as ~ A^Ao <K Ao. 
Therefore, the DQC1 measurements with standard deviation 
Ax refines our knowledge about the interaction strength 9 
from a large standard deviation Ao to a much smaller one 
~ AxAo. By iterating this procedure, the standard deviation Ag 
would decrease successively as Ao — * AxAq — > A x Ao 
With the aid of the adaptive Bayesian inference, Boixo and 
Somma 11 performed a quantitative analysis about this iter- 
ation and concluded that the QML Ag = 0(1 /T) could be 
achieved for an arbitrary desired standard deviation, where 
T — 2 t is the total duration of the estimation process. 

In the next subsection, we discuss the problems of DQC1 
model when it is directly applied to estimate the hyperfine in- 
teraction in the NV center. Before that, we mention a use- 
ful extension (which can be readily verified) of this model: 
the analytical expressions for the quantity estimated by the 
measurement [e.g., (X) = cos(9t) for the DQC1 model and 
(Z) = cos(At) for our protocol, see Eq. ((5)] remains valid for a 
more general initial state p tar = 1 /2+q z o~ z /2 of the target qubit 
with an arbitrary polarization q z . This fact is especially impor- 
tant for estimating the hyperfine interaction in the NV center 
since in this case, initializing the control qubit (the electron 
spin in the NV center) will partially polarize the target qubit 
(the 15 N nuclear spin in the NV center)^ 
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FIG. 2. Energy level diagram of the ground state of an NV center 
in diamond. The hyperfine energy splitting at B = 0.2 T is sketched 
within the \m s = 0) and \m s = -1) manifold. D = 2.87 GHz is the 
zero field splitting of the electron spin, and A is the longitudinal hy- 
perfine interaction to be estimated. The nuclear Zeeman splitting is 
omitted in the diagram. 



B. Direct application of DQC1 parameter estimation to NV 
center 

We consider a negatively charged NV center in diamond 
consisting of a substitutional 15 N atom and a neighboring car- 
bon vacancy. Its electronic ground state is a two-electron spin 
triplet described by a spin-1 operator S, with a zero-field split- 
ting D « 2.87 GHz (described by the term DS 7 :) between 
the \m s = 0) state and the \m s = ±1) states. Under an exter- 
nal magnetic field B along the N-V axis (defined as the z di- 
rection), the Zeeman term g e fisBS z with g e = 2.0023 shifts 
the state \m s = +1) away from the other two states under a 
moderate magnetic field B ~ 0.2 T (see Fig. Thus we 
identify |0) = \m s = 0) and |1) = \m s = -1) as the two states 
of the control qubit of the DQC 1 model and use X, Y, Z as 
the three Pauli matrices for this qubit. The electron spin S 
is coupled to the neighboring 15 N nuclear spin- 1/2 I (with 
the two-fold degeneracy lifted by the Zeeman term gNfi^/BL, 
where g^ = -0.5664 17 ) through the hyperfine interaction 
AS Z I Z + (A ± /2)(S + L + §J + ), where A * 3.03 MHz and 
A j. ~ 3.65 MHz. 17 We regard this nuclear spin-1/2 as the 
mixed-state target qubit of the DQC1 model and use & x , <x v , er- 
as the three Pauli matrices 2I X , 2I y , 2L for this qubit. The diag- 
onal part AS J z of the hyperfine interaction makes the nuclear 
(electron) spin splitting energy dependent on the state of the 
electron (the nucleus). Thus AS J, plays the central role in 
coherent control and readout of the electron and nuclear spin 
states. The hyperfine interaction strength A is the parameter 
to be estimated. 

In the two-qubit subspace, the Hamiltonian H — Hq + H m i x 
consists of the diagonal part 

A) = \gN^ N B(T z + |1><1|<8 (D' - -Afr z ) 
and the off-diagonal part 

£mix = (A J V2)(|04> <1. Tl + |1. T> (0, 
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FIG. 3. (a) Direct application of the DQC1 model to estimate the 
hyperfine interaction strength A in NV center, (b) Combination of 
spin echo and the DQC1 model. 

The diagonal part Hq accounts for the free nuclear spin pre- 
cession with angular frequency g N /i N B, the free electron spin 
precession with angular frequency D' = D - g e [i B B, and the 
projection |1) (1| ® (-Ait z /2) of the diagonal hyperfine inter- 
action AS z l z in the two-qubit subspace. The off-diagonal part 
Hmix is the projection of the off-diagonal hyperfine interaction 
(A ± /2)(§ +1- + S -1+) in the two-qubit subspace. The diago- 
nal hyperfine interaction term |1)(1| ® (-Ach,/2) in Hq cor- 
responds to 77dqci in Eq. ([]]) with 9 <-> (-A/2). It makes 
the precession frequency a> e of the electron spin dependent 
on the hyperfine interaction strength A and the nuclear spin 
state: u> e ^ — D' - A/2 for the nuclear spin state being ||) and 
u> e i — D' + A/2 for the nuclear spin state being |X). There- 
fore, following the procedure in Fig.Q] the interaction strength 
A is encoded as a phase of the electron spin and subsequently 
extracted by estimating (X). 

As schematically shown in Fig. [3 a), the electron spin needs 
to be prepared in the superposition (|0) + |1))/ y2. This can 
be achieved by optical pumping 18 followed by a coherent ro- 
tation. However, this preparation process inevitably influ- 
ences the nuclear spin and changes its state from the unpo- 
larized thermal equilibrium state p" q = 1/2 to a state p[j eq = 
I/2 + q z & z /2 with a finite polarization q z — Ti[p"& z ]^- Then 
the two qubits evolve under the Hamiltonian H for a duration 
r, followed by a measurement of (X). Below we calculate (X) 
without H m j x and then taking it into account by perturbation 
theory. 

Without H m \ x , the two qubits are driven by Hq, which has 
four eigenstates |0, f) , |0, 1) , |1, T) , |1 ; D- The physics is sim- 
ilar to the DQC1 model described in the previous subsection: 

• If the nuclear spin is in the spin-up state |f), then Hq 
drives the precession of the electron spin qubit with an- 
gular frequency io e ^ and the repeated measurements of 
X estimate (X)f = cos(« e f r). 

• If the nuclear spin is in the spin-down state |X), then 
Hq drives the precession of the electron spin qubit with 
angular frequency a> e ^ and the repeated measurements 
of X estimate (X)^ - cos(w c jt). 
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• Now the nuclear spin is in an incoherent mixture of |f) 
[with weight (1 +q z )/2] and ||) [with weight (1 -q z )/2]. 
Then the repeated measurements of X estimate the 
weighted average of (X)i and 

1 + g, ~ 1 — a, j, 
(X) = -^<X> T + -^(Di. (3) 

Then we consider the complications caused by the off- 
diagonal part H m i x . To reduce its detrimental effect on the 
parameter estimation, we consider a suitable magnetic field 
strength (e.g., B = 0.2 T, as indicated in Fig. [2]and used in our 
estimation, see Sec. IIV CI ) so that \D'\ » |AJ. In this case, we 
can use perturbation theory to treat ft m i x , which modifies the 
eigenstates and eigenenergies of the two-qubit Hamiltonian 
H = Hq + H m { x : 

1. ffmix changes the eigenstates of ft from [|0, T), |0, 4), 
|1, T>, IU>] to [|0, T>, |0TT>, I U>, 1 1, 1)1 where 

|6Ti> = [1 - 0(/7 2 )] |0,4) + O(/7)|l,T), 

ii?r> = [i-o( J 7 2 )]ii,t)+o(7 7 )ioa>, 

and T] = AJ(D' + gN^NB - A/2) ~ 10~ 3 for B = 0.2 T. 
In other words, ft m \ x introduces new O(rj) components 
into the eigenstates. It can be readily verified that this 
changes (X) by Oirj 2 ). 

2. H m \ x changes the eigenenergy of |0, |) (|1, T)) by a small 
amount —5 (+6), where 6 = r\A^I2 + 0(r/ 2 A ± ). This in 
turn changes the precession frequencies of the electron 
spin from a> e ^ to (b e ^ = aj ejJ + 6 (jj =T, I). Therefore, 
the average value (X) is obtained from Eq. Q by renor- 
malizing co e41 with cb e4l (yU =f, 1). 

Collecting both corrections discussed above, we obtain 

A 

(X) = cos[(D' + 6)t] cos(-t) (4) 

+ q z sin[(D' + 6)t] sin(-r) + Otf). 

It contains not only A but also undesired parameters such as 
D' (free electron spin precession frequency), 5 (energy shift 
by //mix), and q z (partial nuclear spin polarization). For an ac- 
curate estimation of A, it is desirable to eliminate these unde- 
sired parameters from (X) by modifying the DQC1 protocol. 



which consists of an electron spin n rotation R e v (7r) - e mY ^ 2 = 
-it and a nuclear spin n rotation R"(tt) = e~ OTO >/2 _ 

sandwiched by the free evolution e~' Hr . This composite evo- 
lution contains a spin echo (the part inside the parenthesis) for 
the electron and the nucleus, which eliminates the free pre- 
cession of the electron spin and the nuclear spin. To analyze 
t/com in more detail, we first ignore the off-diagonal part H m [ x 
and then take it into account by perturbation theory. 

Without H m [ x , the Hamiltonian H' = & v YHYfr y commutes 
with H. Thus U com reduces to 

#com = R e y (7T)R';,(n)e- i( "' + " )T = R^n)R';,(n)e- iAT ^ /2 e- ifl ^ T , 

where H ec ho = |1) A& z ) corresponds to //dqci in Eq. (fl} 

with0<-> -A. The operation R"(7i)e^' AT<T ^ 2 on the nuclear spin 
alone can be dropped since it does not influence our measure- 
ment on the electron spin. Therefore, the composite evolution 
becomes ufj m = R1,(7T)e~' HahoT , in which all the undesired pa- 
rameters have been eliminated. 

In the presence of H m \ x , H' consists of the diagonal part 

H' = -jgNMNBfr z + l -Afr z - |1)<1| ®{D' + l -Aa z ) 

and the off-diagonal part H m [ x . Similar to the two-step anal- 
ysis leading to Eq. (0|, H m i x modifies the eigenstates and 
eigenenergies of H — Hq + // m ; x and H' = H'^ + // m ; x : 

1. H m i x introduces new O(rf) components into the eigen- 
states of H and H' . This changes {X) by 0(rj 2 ). 

2. For the Hamiltonian ft, the presence of H m i x changes 
the eigenenergy of |0, 1) (|1,T» by -6 (+6). For the 
Hamiltonian ft', the presence of ft m \ x changes the 
eigenenergy of |0, |) (| 1 , t)) by +6 (-<5). In other words, 
the opposite energy shifts for ft and ft' induced by ft m \ x 
cancel each other in the evolution £/ com . 

For (X), the composite evolution including both corrections 
discussed above is equivalent to 

U com = R^n)e- i(, ^ T + Oirj 2 ), 

i.e., the spin echo eliminates all the named undesired parame- 
ters and the effective evolution U com for the NV center recov- 
ers the DQC1 evolution e _,H DQci T U p to a trivial electron spin n 
rotation R\(n). 



III. ELIMINATING UNDESIRED PARAMETERS BY SPIN 
ECHO 

To remove the dependence on the undesired parameters in 
(X), we combine the DQC1 model with the spin-echo tech- 
nique by replacing the free evolution e~' Hr with the composite 
evolution [see Fig. [2b)] 

U com = e-'^R^R^e-^ = R^R^n) ( e -^M^,)r e -i#rJ ; 



IV. QUANTUM MEASUREMENT PROTOCOL OF 
HYPERFINE INTERACTION 

In this section, first we give the quantum circuit for a sin- 
gle estimation of the hyperfine interaction strength A in the 
NV center. Second, we describe in detail the procedure of 
the entire estimation protocol: the successive adaptation of the 
quantum circuit for dramatically reduced standard deviation 
by combining our prior knowledge with the outcomes of the 
previous measurements through adaptive Bayesian inference. 
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FIG. 4. Quantum circuit for a single estimation of the hyperfine in- 
teraction strength A in the NV center. 



Third, we demonstrate that this protocol approaches the QML 
Aqml = 0(1 /T) for the ideal case. Finally, we include the es- 
sential errors (the nuclear spin rotation error and the electron 
spin decoherence) and show that our protocol still exceeds the 
SQL. 



A. Quantum estimation circuit 

Fig-E]gives the sequence of quantum operations for a single 
estimation of the hyperfine interaction strength A in the NV 
center: 

1 . The electron spin is prepared into the pure state |0) by 
optical pumping. 18 A subsequent n/2 rotation R e x (nl2) 
initializes the electron spin into the superposition |+) = 
(|0> + |1»/V2. The nuclear spin is a partially polar- 
ized state p„ = 1/2 + q z & z J2. This initial density matrix 
Anitiai = ® Pn coincides with the initial density 
matrix Pdqci = 1+) (+1 ®Aar of the DQC1 model, where 
the target qubit state p tal - also has an arbitrary polariza- 
tion, as discussed at the end of Sec. Ill Al 

2. The two qubits experience a composite evolution 
(within the dashed box in Fig. @}, which consists of a 
free evolution e~' Hr , a controlled nuclear spin n rotation 
R n y (n) = ® (-ifry) + |0> <0|, an electron spin ro- 
tation Ri(n), another controlled nuclear spin n rotation 
R"(n), and another free evolution e~' Hr . The equality 
R"/n)R*(7T)R'l,(7r) = R^(n)k e y (n) shows that this compos- 
ite evolution coincides with t/ com in SecUTA 

3. A 7r/2 rotation Rl(n/2) is applied to the electron 
spin, followed by a measurement of Z through optical 
methods. 4,19 This measurement estimates 

(Z) = TrZ^(^/2)f> com p lllitia if/] om [^(7r/2)] t 



Since the evolution e lH ^° T 
same form as the DQC1 model, the average value is 



e -iH DQC ,T| fl ^_ A has t he 



<Z> = cos(Ar) + O(rf). 



(5) 



The electron spin rotation R e x (nl2) [R e y (n)] in the circuit 
is achieved by a tt/2 pulse (n pulse) with the central fre- 
quency \D'\ and the bandwidth » A/2, so that both transitions 
|0,T> <-> |1,T> and |04> <-> |1,|> are equally excited. The 



controlled nuclear spin rotation R"(n) is achieved by a n pulse 
centered at the resonant frequency A-gNfi^B -6 of the transi- 
tion 1 1 , t) — > |1» D- The duration t of the free evolution can be 
chosen in the experiment as t > 1/A ~ 0.1 fis. The electron 
spin rotation occurs within a few nanoseconds and hence can 
be regarded as instantaneous^Si^l However, the controlled nu- 
clear spin 7r rotation takes t„ ~ a few microseconds, compara- 
ble to the free evolution time t. Detailed analysis in appendix 
A shows that incorporation of t„ amounts to replacing the free 
evolution time r in Eq. Q by the sum (r + r„). For brevity, 
we use t to denote (t + r„) from now on. 

In arriving at Eq. (0, we have assumed that all the gate op- 
erations in the circuit and the measurements of Z are free of 
errors. In a realistic experiment, the most basic errors include 
the deviation of the nuclear spin rotation angle from n in the 
controlled n rotation R"(x) and the finite electron spin coher- 
ence time T e 2 : 

• Nuclear spin rotation error. The two controlled nu- 
clear spin n rotations R"(n) in the quantum estimation 
circuit (Fig. HJ are subjected to random errors, which 
may come from our limited prior knowledge (which be- 
comes more and more precise after each successive es- 
timation step) about the interaction strength A or other 
experimental sources. For the actual rotation angle 
(n + 2e) differing from n by an error 2e, the actual con- 
trolled rotation R"(tt, e) = R"(n) + 6"(n) differs from the 
ideal one R"(n) by 

~6 n y (n)=\l)(l\®(-e+i^& y ) + 0(<?). 

For the first controlled rotation being R"(^, £ a ) an d the 
second controlled rotation being R"(n, the actual 
quantity estimated by the quantum circuit M(t) is 



<Z e ) 



1 



(4) + <4A 



cos(Ar) 



+ (e a e b ) + (e a )0(i]) + (e h )0(i]) + 0(if), 



The first source of error is our ignorance about A. 
In the £-th estimation step, our limited prior knowl- 
edge about A (as quantified by the standard deviation 
Ajt-i, see Sec. IIV Bl and hence the resonant frequency 
A - gN^NB - 6 of the transition |1, f ) — > \l,i) makes 
it impossible to construct an exact n pulse for this tran- 
sition. The typical detuning for this transition is Ajt_i. 
The typical rotation angle deviates from the ideal value 
n by an amount nt^ k _ l /(2Q. 2 ) ~ 10~ 3 , the same order of 
magnitude as 0(rj), for the Rabi frequency Q = 500 kHz 
used in our estimation. Thus every term in the second 
line of the above equation has the same order of ~ 10~ 6 , 
which allows us to replace the second line by 0(rf). For 
€ a and eh being independent, we obtain 



(Ze) = (1 



e 2 ) cos(At) + 0(?] 2 ), 



where s 2 = (e 2 ) = (e?). For other experimental sources, 
the errors are typically random with (e a ) = (e/,) = 
(£a£b) = 0, so that the above equation still holds. 



6 



Electron spin decoherence. The electron spin in the 
NV center is subjected to decoherence by the surround- 
ing 13 C nuclear spin bath. The coherence time of the 
electron spin in the ground state is Tt ~ 350 /^s under 
the natural abundance of the 12 C isotope (98.8%), and 
it is extended to 1.8 ms under the ultrapure 12 C abun- 
dance (99.7%) at room temperature^^ By incorpo- 
rating the electron spin relaxation (with the relaxation 
time 23 T* =5.9 ms) and decoherence in the Lindblad 
form, it is straightforward to show that the quantity es- 
timated by the quantum circuit is no longer Eq. (0 but 
instead 



(Z d ) = e- 2T/T "i cos (At) + 0(t] 2 ). 



In summary, in the presence of errors, the quantity esti- 
mated by the quantum circuit in Fig. [4]is given by 



<Z) = e(r)cos(Ar) + (9(77 2 ), 



(6) 



where Q(j) — 1 - e 2 for the nuclear spin rotation error of 
magnitude s and Q(j) = e" 2T ^- for a finite electron spin co- 
herence time r|. In our estimation, we use B = 0.2 T so that 
the correction for the hyperfine interaction Oirj 2 ) ~ 10~ 6 . 



B. Estimation procedure 

We use M(t) to denote the quantum estimation circuit in 
Fig. |4] whose total duration is 2t. A single run of the circuit 
M(t) returns two outcomes: +1 for the electron spin in the 
state |0) or -1 for the electron spin in the state |1), with cor- 
responding probabilities p ± \ — [1 ± (Z)]/2. An estimator of 
the average value (Z) [Eq. ©J is obtained by averaging over 
the outcomes of repeated running of the circuit. For example, 
averaging over N measurements produces Z, a single estima- 
tor of (Z). By the central limit theorem, for relatively large N 
(e.g., > 100), this estimator obeys the Gaussian distribution 
N({Z), £) centered at (Z) with a standard deviation C, — 1 / "0V. 
Alternatively, we can also say that the average value (Z) obeys 
the Gaussian distribution N{Z, (), which actually means that 
the difference (Z) - Z obeys the Gaussian distribution A/(0, (). 

The estimation begins with a prior knowledge of the hy- 
perfine interaction strength A. It is quantified by a Gaussian 
distribution AAAo, A ( >) centered at Ao with a relatively large 
standard deviation Aq, which quantifies our ignorance about 
A. This prior knowledge tells us, with a 95% confidence, that 
A lies within the interval [A - 1.96A ,A + 1.96A ]. From 
the prior knowledge A/(Ao,Ao), we construct the quantum 
circuit M(j\) for the first estimation, which provides a new 
knowledge about A, as quantified by a Gaussian distribution 
N{A\, Ai). Through the Bayesian inference, this new knowl- 
edge is combined with the prior knowledge to produce an up- 
dated knowledge about A, quantified by a Gaussian distribu- 
tion N{A\,A\) with a smaller standard deviation Ai < Ao. 
Therefore, the first estimation step refines our knowledge 
about A from N(A , A ) to AAAi, Ai) (with Ai < A ), which 
in turn serves as the prior knowledge of the next estimation 



step. By iterating this procedure, the standard deviation of 
the Gaussian distribution quantifying our ignorance about A 
would decrease successively as Ao > Ai > A2 > ■ ■ ■ . The 
iteration is stopped at the K-th step when the desired standard 
deviation Ajesire is achieved: A^ < Ad es ire- Below, we describe 
the above estimation procedures in more detail. 



1. Gaining knowledge about A from measurements 

In the k-th estimation step (k = 1 , 2, ■ • • ), the prior knowl- 
edge about the hyperfine interaction strength A is quantified 
by the Gaussian distribution N(A k -i, A k -[). Suppose that T k 
has been properly chosen (to be discussed shortly). By run- 
ning the circuit M(j k ) for a relatively large number N k (> 100) 
of times, we obtain an estimator Z k of {Z) k = Q(r k ) cos(Ar,t) + 
0(t] 2 ) with a standard deviation & = 1/ V^Vjt- This knowledge 
tells us that {Z) k obeys the Gaussian distribution M(Z k ,( k ). 
We need to convert this distribution of {Z) k to a distribution 
of A. For a general T k , the relation between (Z) k and A is 
nonlinear and the conversion from (Z) k to A results in a non- 
Gaussian distribution of A, with a characteristic width 



fi* 



fit 



\d{Z) k /3A\ Q(T k )r k \ sin(AT jt )| 
Now we determine according to two requirements: 

1. The distribution of A should be Gaussian (i.e., the 
relation between {Z)u and A should be linear), so 
that analytical results can be obtained. Based on our 
prior knowledge N(Ak-\, An) about A, the conditions 
Ak-\Tk — tt/2 + 2jt x integer and Ak-iTk ^ 1 enable the 
Taylor expansion {Z) k = (A k ^ -A)Q(r k )T k + S k + 0{rf) 
with** « Q(T k )(A k - lTk ) 3 /6. For 6 k ,0(?] 2 ) « ( k ,\(Z) k \, 
the correction terms 6 k + 0{rj 2 ) can be safely dropped, 
so that the relation between (Z) k and A becomes linear 
and the distribution of A becomes Gaussian N(A k , A k ) 
with 



Ak 



Q(T*)T* 

a 1 



Q(T k )T k Q{T k )T k Wk 



(7a) 
(7b) 



The distribution N(A k ,A k ) of A tells us, with a 95% 
confidence, that A lies in the interval [A k - l.96A k ,A k + 
1.96 A,]. 

2. For maximal precision of the estimation, the standard 
deviation A* should be minimized, i.e., Q(T k )T k should 
be maximized. 

Eq. dTbl shows that the standard deviation At of the mea- 
surement of A is equal to the standard deviation fit = 1/ VA^ 
of the measurement of (Z) k divided by Q(j k )T k : 

• For Q(T k ) = 1 (i.e., no errors), the standard deviation 
A k is reduced upon the increase of T k , which can be 
interpreted as a repetition of the circuit operations (as 
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enclosed in the dashed box in Fig. |4|i before the mea- 
surement is made. This is equivalent to a multiround 
protocol suggested by Giovannetti et al? . Therefore, 
the dependence A* oc 1 /tj implies the QML. 

• The standard deviation A\ is reduced upon the increase 
of Nk- The dependence A* oc 1 / -\/Nk implies the SQL. 

In summary, for optimal performance, we should first 
choose (or equivalently Nk) subjected to the constraint 



0(if) « Ck « 1 



(8) 



and then choose Tk to maximize Q(Tk)Tk, subjected to the con- 
straints 



At-iTjt = - + 2m k n, 



i^k-lTkY 



Q(T k )Ak-iT k » O(if), 



(9a) 

(9b) 
(9c) 



where m k e Z and 0(rf) ~ 1(T 6 for B = 0.2 T. The con- 
straint (k <Z 1 ensures the validity of our Gaussian distribu- 
tion assumption for (Z)k, while other constraints ensure the 
validity of the formula (Z)j ~ (A*_i - A)Q(jk)Tk- The er- 
ror of the linear expansion can be dropped if 5k <s which 
gives Eq. ( |9bl with Q(r^) < 1. Eq. d9cT > denotes the condi- 
tion to drop 0(tj 2 ) in (Z)k- Note that the constraints [Eqs. (0] 
on Tk have no solution under certain conditions, e.g., when 
Q(Tk) ^ Ofjfy/fj^k) '■ Therefore, for more flexible choice 
of Tk, the standard deviation ft of the measurement of (Z)t 
should not be too small. 



For very accurate measurement compared with the prior 
knowledge, i.e., A* <s Aj_i, Eqs. fllOal ) and d 1 Obb reduce to 
Ak ~ Ak and Ak ~ Ak, suggesting that the updated knowledge 
is dominated by the measurement. By contrast, for inaccurate 
measurement » At-i, the updated knowledge At w Ak-i 
and Ak ~ At-i is dominated by the prior knowledge. 



C. Ideal case: approaching quantum metrology limit 

In this subsection, we demonstrates the QML scaling of 
our estimation protocol in the ideal case, i.e., in the absence 
of any errors (e.g., operation errors, relaxation, and decoher- 
ence). For simplicity, we assume that in each estimation step, 
we run the quantum circuit for the same number of times Ni = 
N 2 = ■ ■ ■ = N, corresponding to fi = & = • • • = f = 1/ V^V- 

Up to the K-th estimation step, the total duration of the our 
estimation process (identified as the total amount of resources 
spent) is 

K 

R K = N^2T k = NTf. 

t=l 

To see the scaling of the precision 1/Ai with respect to Rk , 
we take the first estimation step as a reference. Further, we 
take Ao = °o to exclude the contribution from the prior knowl- 
edge N(Ao, Ao), so that all our knowledge about A comes from 
the resources Rk spent in our protocol. Then, the QML limit 
Ajc,qml is defined by A^qml/A! = l/(R K /R{), while t he SQL 
limit A/t,sql is defined by Ajf,sQ L /Ai = 1/ \JR K /Ri- Using 
Rj = 2Nt x and Ai = l/(r, y/N), we obtain 



2. Combining new knowledge with prior knowledge 

In the previous subsection, we have spent Nk runs of the cir- 
cuit M(jk) to obtain the new knowledge N(Ak, Ak) about A. To 
make use of the resources spent in obtaining the prior knowl- 
edge N{Ak-\, Ak-\), we use the Bayesian inference, which 
combines our new knowledge N(Ak, Ak) with the prior knowl- 
edge M(Ak-\, Ak-i). It gives an updated Gaussian distribution 
N(Ak, Ak) centered at 



Ak-JA 2 ^ +A k /A 2 k 



(10a) 



(which is a weighted average of Ak-\ with weight 1/A?_, and 
Ak with weight 1/A?) with a standard deviation A* determined 
by 



a; 



A 2 

"t-1 



t2" 



(10b) 



This updated knowledge N(Ak,Ak) tells us, with a 95% con- 
fidence, that A lies in the refined interval [Ak - 1.96 At, At + 
1.96At]. The inequalities At < An and At < reveal that 
the combination of N{Ak-\, A^) and A/(At,At) gives us a 
more precise knowledge about A. 



1 



A 2 

A",QML 



= N 



A 2 



I 

t=i 

K 

NTi ^ T k , 



Tk 

s k=\ ' 
K 



(ID 

(12) 



k=\ 



1 K 
A^= N ^' 

k=\ 



First, we compare Ak with the QML limit A^qml and the 
SQL limit A^sql and discuss the condition for approaching 
the QML: 

1. The inequality Ak > A^qml can be readily verified. 
This manifests the QML precision 1/A^.q ML as the up- 
per precision bound. To achieve the QML, {ft} should 
satisfy tk » tk-i » • •• » ti, so that the total amount 
of resources is dominated by the final estimation step 
and hence Ak ~ A^qml ~ 1 / (tk yN). This condi- 
tion is equivalent to a dramatic reduction of the stan- 
dard deviation of the measurement for each successive 
estimation step: A^ <k Ak-i <k ■ • • <c A t . This ensures 
that in each estimation step (say, the k-th step), the stan- 
dard deviation of the estimation, At w At w l/(rt yfR), 
is dominated by the standard deviation A* of the mea- 
surement instead of the standard deviation At-i * At_i 
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FIG. 5. Comparision of the standard deviation A K of our protocol 
with the QML limit A^qml (solid line) and the SQL limit A^sql 
(dashed line). How to choose r* is explained in the main text. 



of the prior knowledge [cf. Eq. ( UObH . The condition 
Tk » tk-\ s> ■ • • » t\ is also equivalent to 

Afc.iTjfc^, (13) 

since A H r ( w i-rt = (r k /T k -i)C 

2. For ti = T2 = • • • = Tfc, the precision = NKt 2 co- 
incides with the SQL precision 1/A^, S q L since in this 
case our protocol reduces to simple repetition of the 
same quantum circuit M(t\). 

Then we give the best choice {r^ dea1 } satisfying the QML 
condition Eq. ( fT3l ) for the ideal case according to the descrip- 
tion in Sec. IIV Bl We choose {r^ dea1 } by taking the largest 
such that Aj_iTjt ~ c at every step, where c is a constant satis- 
fying c » f and c 3 <K 6£. Then {T^ dea1 } automatically satisfies 
the QML condition Eq. ( fT3l and the linear expansion condi- 
tion Eq. ( f9bb . From A^-iT* ~ c, we have Ak w A* « (£/c) k Ao, 
i.e., the standard deviation A* ~ A* is dramatically reduced by 
each successive estimation step. We also have r^ 63 ' w (c/%) k TQ 
(with To defined through Ao = (/tq), i.e., an exponential in- 
crease of T^ deal with k. Note that, for B = 0.2 T, we have 
0(t] 2 ) ~ 10- 6 . Therefore f can be as small as ~ 10 . 

Finally we provide a numerical simulation for the estima- 
tion process. The parameters for the simulation are A = 
3.06 MHz, B = 0.2 T, N = 1000, corresponding to f « 0.03. 
We take c = 0.2, which satisfies c » g and c 3 «: The 
prior knowledge is Ao = 3.03 MHz with a standard devia- 
tion Ao = 0.03 MHz, which has been reported by a previous 
experiment. 17 Each controlled nuclear spin n rotation uses a 
l-/is square pulse with the Rabi frequency Q = 500 kHz. The 
electron spin rotations are regarded as instantaneous, as men- 
tioned at the end of Sec. IIV A| In Fig. [5] the proximity of A^ 
(circles) to A^qml (solid line) confirms the QML scaling of 
the estimation. 



D. Realistic case: surpassing standard quantum limit 

In this subsection, we take into account the nuclear spin 
rotation error and electron spin decoherence and discuss the 
optimal choice of {t^} and the resulting precision 

^=/vf>(T*)T t ] 2 

k=\ 

of the estimation, derived from Eq. (fTbl and (llObl i: 

• Nuclear spin rotation error Q(t) = 1 - s 2 = Q. This 
error is equivalent to an increase of £ to \ = £/ Q. Then 
QML condition Eq. ( fT3l becomes Afc_iT* » if. For a 
general Q that is not too small (i.e., 1 > Q » £), the 
conclusion in the ideal case remains valid with £ — > <f, 
i.e., {rit} is chosen as ~ (cICfijalQ), where c is a 
constant subjected to c » Az and c 3 <s. 6Az- In the 
simulation, we consider a typical error s — 0. 1 (corre- 
sponding to ~ 3% error in the rotation angle). Then we 
have Q ~ 1, and this allows us to set c = 0.2, the same 
value with the ideal case. As a result, we can choose 
Tk ~ i"J, deal and Ak is nearly the same as the ideal case. 
Therefore, the QML scaling is preserved for the realistic 
nuclear spin rotation error, as confirmed by the nearly 
complete coincidence between Ak (stars) and A^qml 
(solid line) in Fig. [5] 

• Electron spin decohence Q(t) = e~ 2r t T i . According to 
Sec. IIV B| we should choose to maximize Q{Tk)Tk, 
subjected to the constraints in Eqs. (0. We use Ak-\Tk ~ 
c = 0.2 in the simulation. In the presence of the electron 
spin decoherence, Q(t) decreases as t increases. Thus 
the QML condition c » (IQ(T k ) is no longer valid at 
some point. This is why A* starts to deviate from the 
QML line at k = 3 in Fig. [5] Note that the estimation of 
k — 3 still surpasses the SQL. The maximum of Q(j)t 
occurs at r = T e 2 /2, meaning that the standard deviation 
Ak of the quantum circuit M(jk) is the smallest when 
T k ~ 7*1/2. Further increase of Tk makes the precision 
of M(jk) worse. Once Tk reaches Tk ~ TZ/2 at k 

k, . the estimation for k > k c is performed with Tk = 
Tk c . Therefore, for K — k c + K, further estimation steps 
beyond k c (i.e., k = k c + 1, • • • ,k c + K) increases the 
precision 1 / A 2 K by the SQL trend: 

- -i- * NK(T e 2 /2) 2 . 
A 2 . A, 

k c +K K 

For r| = 350 /is, we have k c = 4. Fig.|5]shows that Ak 
surpasses the SQL for K < 4, while it decreases parallel 
to the SQL for K > 4. 



V. CONCLUSIONS 

We have proposed an efficient quantum measurement proto- 
col to estimate the hyperfine interaction between the electron 
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spin and the N nuclear spin in the NV center. The essential 
idea of our protocol is the combination of the DQC1 parame- 
ter estimation 11 with the spin-echo technique. The spin echo 
eliminates the independent dynamics of the electron spin and 
the nuclear spin in the DQC1 model, but keeps the dynamics 
due to their interactions, whose strength is to be estimated. 
This protocol does not require the preparation of the nuclear 
spin state. We quantify the resources R as the total duration 
2 t of the estimation process. In the absence of any errors, the 
precision 1/A 2 (with A being the standard deviation) of the 
estimation approaches the quantum metrology limit (QML) 
1/Aq ML = 0(R 2 ). This QML scaling is robust against the typ- 
ical nuclear spin rotation error in realistic experimental con- 
ditions. In the presence of electron spin decoherence, the pre- 
cision 1/A 2 keeps its QML scaling when r <K TZ/2. Once r 
becomes close to T% further estimation steps increase the pre- 
cision 1/A 2 according to the scaling 1/A 2 ,^ = 0{R) of the 
standard quantum limit (SQL). Due to the QML scaling in 
the initial stage, the overall precision still surpasses the SQL. 
We expect that this method can be applied to other solid state 
systems such as quantum dots or cold atoms to measure the 
interaction between two spins. 
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Appendix A: Accounting for finite duration of controlled 
nuclear spin rotation 

In this section, we assume that each of the two controlled 
nuclear spin n rotation in the quantum protocol (Fig. |4) is 
driven by a square n pulse with a duration r„ ~ 1 fis and 
prove that inclusion of this finite duration amounts to a trivial 



renormalization r — > t + t„ in Eq. ©. 

In Fig. |H the initial state pi n j t i a i = |+) (+| ® p n is prepared 
at t = -t — t„. The first free evolution e~ lHr occurs during 
t € [— t - t„, -t„], followed by a controlled nuclear spin n 
rotation during t e [— r„, 0]. A fast electron spin n rotation 
is applied at t — 0, another controlled nuclear spin n rotation 
during t e [0, t„], and another free evolution e~' Hr during t e 
[t„,t + t„]. 

First we calculate the evolution operator driven by a square 
ji pulse applied during t € [fi.fc], with a central frequency 
a> = A - gNHfjB-5 (where +6 is the energy correction to |1,T) 
by the off-diagonal hyperfine interaction) resonant with the 
transition |1,T) — > ILD- During this pulse, the Hamiltonian 
H(t) = H + V(t) of the electron-nuclear spin qubits acquires 
an additional term 

V(t) = l ^-{e~ iM 1 1 4) < 1 , Tl - e iM 1 1, T) <1, 41), 
with a constant Rabi frequency Q. R = n/fa - fi). With the 
aid of the interaction picture ^(f)} = e ,H '| v F(f)), the evolu- 
tion operator Uv(t2,h) during t e {t\,t%\ can be calculated 
straightforwardly as U v (t 2 ,t\) = e~ i & h e~ i6l<h ~ >l ' ) eP th , where 
Hi(t) = e'"'V(t)e ~ iftt . Similar to the discussions in Sec. Ill Bl 
we have Hi(t) = (Cl R /2)\l)(l\®& y + 0(Q. R 7]), where rj ~ 10~ 3 
for the external magnetic field B = 0.2 T used in our esti- 
mation. Therefore, the evolution e~ lH ^ }1 ~ h) a R"(t[) coincides 
with the instantaneous controlled rotation and hence 

U v (t z ,h) = e"*^)/' 1 . 

With the aid of this result, it can be readily checked that the 
evolution operator for the composite evolution (as enclosed by 
the dashed box) in Fig. |4]is equal to t/ com | T ^( T+Tn ). Therefore, 
inclusion of the finite duration t„ of the controlled nuclear 
spin rotation amounts to replacing t with (t + t„) in Eq. ([5]). 
Note that the nuclear spin relaxation time and decoherence 
time > 1 ms are much longer than r„ ~ 1 jus and hence have 
negligible influence on this result. 2 ^ 
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